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(N) ABSTRACT 
CO 

We present 12 CO(2-l) observations of the submillimetre galaxy ALESS65.1 per- 
f--^ formed with the Australia Telescope Compact Array at 42.3 GHz. A previous ALMA 

study of submillimetre galaxies in the Extended Chandra Deep Field South detected 
j— 1 [CII] 157.74 /im emission from this galaxy at a redshift of z — AAA. No 12 CO(2-l) emis- 

L* sion was detected but we derive a firm upper limit to the cold gas mass in ALESS65.1 

. £^ of M gas < 1.7 x 10 10 M . The estimated gas depletion timescale is < 50 Myr, which is 

similar to other high redshift SMGs, and consistent with z > 4 SMGs being the likely 
$_i progenitors of massive red-and-dead galaxies at z > 2. The ratio of the [CII], 12 CO and 

far-infrared luminosities implies a strong far-ultraviolet field of Go > 10 3 , as seen in 
Galactic star forming regions or local ULIRGs. The observed £[cii]/-^fir = 2.3 x 10 -3 
is high compared to local ULIRGs and, combined with £[cn]/£co ^ 2700, it is con- 
sistent with ALESS65.1 either having an extended (several kpc) [CII] emitting region 
or lower than solar metallicity. 

Key words: galaxies: evolution, galaxies: formation 



1 INTRODUCTION 

Submillimetre galaxies (SMGs) are ultraluminous, dusty 
starbursting systems with extre me star formatio n rates 
(SFRs) of 100 - 1000 Mq yr" 1 (e.g. |Blain et al.|2002[ ). These 
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SMGs have typical redshifts of z 
[2005 



2-3 (e.g. Chapman et al. 
|Wardlow et~!iri|2l)TT| |Yun et al.||2012| |Smolcic et al. 



[2012|, but an increasing number of higher redshift SMGs are 
being found. The z > 4 SMGs represent candidates for the 
most intense star formation phase of the massive red galaxies 
seen at z > 2 ( |Cimatti et al.||2008| |Marchesini et al.|2010 |. 
This high redshift tail of the SMG distribution was long 
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under-represented in SMG redshift surveys, mostly because 
they lie below the sensitivity limits of the radio interferome- 



ter surveys used to identify SMGs. However, in the last 
years, over a dozen of these sources have been reported ( 


few 
Ca- 


pak et al. |2008[ |2011| Daddi et al. |2009a|b| |Coppin et al. 


2009[|Knudsen et al.|2010||Carilli et al.|2010[ 


2011 


|Riechers 


et al.||2010| |Smolcic et al.||2011 Cox et al.| 


2011 


Combes 


et al.||2012 Walter et al.||2012i, and recent ALMA redshift 



surveys have doubled these numbers (Weiss et al. submitted; 
Vieira et al. submitted). Detailed studies of the star forma- 
tion and gas content in these z > 4 SMGs are still rare, 
but are needed to provide a unique insight into the earliest 
phases of the growth of massive elliptical galaxies. 

The most accessible tracer of cold gas in galaxies is 
12 CO, which has been observed extensively in SMGs (e.g. 
Greve et al.|2005||Tacconi et al.|20"06}|Bothwerl et al.|2013[ ). 
Detections of CO in z > 4 SMGs have shown they are gas 
rich systems with sufficient reservoirs (M gas > 10 10 M ) to 
sustain the extreme star formation rates of ~ 1000 Mq yr _1 
for only short time scales (10s of Myrs) ( jCoppin et al.|2010| 
Riechers et al. 20101, unless the gas is replenished. This is 



consistent with high redshift SMGs being the progenitors of 
the luminous red galaxies seen at z > 2. 

The 2 -P3/2 _ 2 Pi/2 fine structure line of singly ionised 
carbon at 157.74 fj,m (here-after [CII]) has emerged as a 
powerful alternative line for studying the ISM in high red- 
shift sources. It can represent up to 1% of the bolomet- 



ric luminosity of star forming galaxies (e.g. Crawford et al. 
1985 Stacey et al. 19911. This line emission arises mainly 



from the photodissociation regions (PDRs) that form at the 
edges of molecular clouds illuminated by the UV photons 
of young-massive stars, but a significant contribution can 
also come from HII regions and the more diffuse warm in- 



terstellar medium (Madden et al. 1993 Heiles 19941. The 



[CII] line therefore provides an important probe of the gas 
content and star formation processes in a galaxy. 

A recent ALMA Cycle study of 126 submillimetre 
sources located in the LABOCA Extended Chandra Deep 
Field South (LESS, |Weifi et al.|[2009| Karim et al. 2012; 
Hodge et al. 2013 submitted) resulted in the serendipitous 



identification of [CII] line emission from two SMGs (Swin- 
bank et al^|2012| hereafter S12). The average [CII[/far- 
infrared luminosity ratio of these two SMGs is ~ 0.0012 
± 0.0004, roughly ten times higher than that observed in 
local ultraluminous infrared galaxies, which has been in- 
terpreted as evidence that their gas reservoirs are more 
extended (S12). The large extent of SMGs is supported 
by other observations, such as extended radio morpholo- 
gies (e.g. Chapman et al. 2004 Biggs & Ivison 20081, ex- 
tended H-a morphologie s (e.g . Swinbank et al.||2006l, large 



CO(1-0) sizes (Ivison et al 



2010a 



2011 



Thomson et al. 



2012 Hodge et al.|2012 i, lack of silicate absorption in mid 



infrared spectra ( jMenendez-Delmestre et al.|2009 l, and less 
reddened broad-band mid-infrared c olours (|Hainline et al. 
20091). Hi g h [CII] /far-infrared ratios ( jStacey et al.|2010||Ivi- 
son et aL||2010b S12) add to this mounting evidence that 
star-formation in SMGs takes place in a region larger than 
the compact nuclear starbursts of local ULIRGs. 

In this paper we present 12 CO(2-l) observations of one 
of the ALMA detected SMGs, ALESS J033252. 26-273526. 3 
(hereafter ALESS65.1). We adopt the standard A-CDM cos- 



mological parameters of £2m = 0.27, Qa = 0.73, and a Hub- 
ble constant of 71 km s -1 Mpc -1 throughout this paper. 



2 OBSERVATIONS AND RESULTS 

The 12 CO(2-l) line (^ rcst = 230.538 GHz) in ALESS65.1 
(RA(J2000) = 03 32 52.26, Dec(J2000) = -27 35 26.3) (S12) 
was observed over a period of four consecutive nights, 9 - 
12 August 2012, with the Australia Telescope Compact Ar- 
ray ( ATC A) , using the Compact Array Broadband Backend 
(CABB). The array was in the most compact five-antenna 
configuration, H75, which has a maximum baseline of 89m 
and two antennas set along a northern spur. The hybrid 
configuration allows good (u, v) coverage to be obtained for 
integrations less than the full 12 hour synthesis. We centered 
the 7mm IF1 receiver on 42.343 GHz, the expected frequency 
of the 12 CO(2-l) line emission given the [CII] redshift of 
z = 4.4445 derived by S12. The 2GHz wide bandwidth of 
CABB results in a frequency coverage of 41.3 to 43.3 GHz, 
covering 12 CO(2-l) emission between z = 4.32 - 4.58. The 
weather was good to average, with rms atmospheric path 
length variations of 100 to 400 microns throughout the run, 
as measured on the 230m baseline ATCA Seeing Monitor 
(Middelberg et al. 2006). The system temperature was 140 
to 250 K throughout the four nights. Weather and atmo- 
spheric conditions can induce temporal fluctuations across 
the wide CABB band, so, following |Emonts et al. (20111, 
a bandpass calibration scan was acquired at the beginning, 
middle and end of each 8 hour night. Phase and amplitude 
calibration information was acquired with 2 minute scans 
on PKS 0346—279 every 10 minutes and pointing checks 
performed on the same source every hour. Flux density cal- 
ibration was performed on Uranus at the beginning of the 
nights, when it was at an elevation of roughly 55 degrees. 
The uncertainty in the flux density calibration using the 
standard MIRIAD model of Uranus is estimated to be 30% 



(Emonts et al. 20111 



The data were calibrated, mapped and analysed using 



the standard MIRIAD ( |Sault fc Killeen 1999) and KARMA 
( Gooch|1996 1 packages. The synthesized beam from natural 
weighting is 14.0 x 9.9 arcsec. A total of about 20 hours 
on-source integration time was obtained over the 4x8 hour 
nights. ALESS65.1 was not detected in the 42.3 GHz con- 
tinuum map from the full CABB band, which achieves an 
rms noise level of 11 /iJy beam -1 . 

The resultant channel noise in the 1 MHz (7.1 km s -1 ) 
wide spectrum is ~ 0.43 mjy beam -1 , consistent with other 
comparable 7mm ATCA/CABB surveys (e.g. |Coppin et ah] 
2010 Emonts et al.|20lT l. The visibilities were resampled to 
velocity resolutions of 50, 100 and 200 km s -1 and each cube 
was examined for an emission line near the ALMA position. 
A ~4er spike at the position of ALESS65.1 and offset from 
the [CII] emission by ~ 280 km s _1 was examined in detail. 
It was deemed a noise spike because it is only present in one 
channel in all three binned images and is present only in 
the last night of data when each night is imaged separately. 
We flag the spike (20 native channels) in the last night of 
data and recombined all four nights of data to derive the 
final spectra. The spectra at the source position in the 100 
and 200 km s -1 binned cubes (Figure 1) have an rms of 
0.11 and 0.08 mjy beam -1 , respectively. Assuming a line 
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Figure 1. 12 CO(2-l) spectrum of ALESS65.1 binned into 100 
km s" 1 channels and extracted at the ALMA position. The red 
dashed line shows the ALMA [CII] spectrum from S12, binned to 
a similar channel width (120 km s — 1 ) and with the flux density 
divided by 15 for clarity. The blue dotted line shows the spectrum 
binned into 200 km s -1 for maximum sensitivity. We conclude 
that no line was detected to a 3<r limit of 0.24 mjy beam -1 per 
200 km s -1 channel. 



FWHM identical to the [CII] emission line (470 km s" 1 , 
S12) and adopting a 3<r limit, the 12 CO(2-l) line intensity 



of ALESS65.1 is I, 



CO(2-l) 



< 0.11 Jy km s 



3 DISCUSSION 

The observed and derived properties of ALESS65.1 are sum- 
marised in Table 1. To estimate the cold molecular gas con- 
tent in ALESS65.1 we calculate upper limits to the line lu- 
minosity and total cold gas (Kb + He) mass from the CO(2- 
1) flux density limit. Following Solomon and Vanden Bout 
(2005), the line luminosity upper limit is £co(2-i) < 2.2 
x 10 10 K km s" 1 pc 2 . If we assume the gas is thermalised 
(i.e. intrinsic brightness temperature and line luminosities 



are independent of J transition), so L' co , 



Li- 



3(2-1) — ^CO(l-O)! 

and a CO-to-Hb conversion factor a = 0.8 Mq (K km s" 1 
pc 2 ) -1 , which is appropriate for ULIRGs (e.g. Downes & 



Solomon|1998 but see Bothwell et al.|2012 I, the upper limit 
on the total cold gas mass is M gaa < 1.7 x 10 10 Mq. This 
is consistent with the gas mass found in other z > 4 SMGs 



(Schinnerer et al. 


2008 


Daddi et al. 2009a Coppin et al. 


|2010||Walter et al. 


2012 


1 and the typical [CII]/M gas ratio at 



high redshift (S12). 

The total baryonic mass of ALESS65.1 can be derived 
by combining the gas and stellar mass estimates for the sys- 
tem. The stellar mass of the system was estimated from the 
rest-frame absolute H-band magnitude to be M* ~9x 10 10 
M Q (S12), so the gas fraction is modest with M gas /M* < 0.2. 
The total baryonic mass A'/bary = M gas +M* is ~ 9 - llx 10 10 
Mq . This is consistent with the dynamical mass for the sys- 
tem, based on the [CII] linewidths and spatial extent, of 
Md yn sin 2 (j) ~ (3.4 ± 1.8) x 10 10 M (S12). 

ALESS65.1 is detected in the 870^im continuum by 
ALMA and weakly detected in the far-infrared (FIR) by 
Herschel (S12). A fit to the IR SED resulted in a restframe 
IR (8 - 1000 ^m) luminosity L m = (2.0 ± 0.4) x 10 12 L Q 



(S12), which corresponds to a st ar formation rate of ~ 340 
Mq yr _1 using the conversion of Kennicutt (19981. The gas 



depletion timescale r = M gas /SFR < 50 Myr is similar to 
the gas depletion rates of other high redshift SMGs ( Schin- 
nerer et al.||2008 Coppin et al.|[2010 1 . This short timescale 



provides further evidence that z > 4 SMGs have the gas 
consumption timescales necessary to be the progenitors of 
red-and-dead "ellipticals" at z > 2. 

Next we examine the physical conditions of the gas in 
ALESS65.1 using the [CII] detection and 12 CO(2-l) limit. 
The L[ CII ]/L F ir versus Lco(i-o)/£fir diagram is a pow- 
erful diagnostic as these two ratios are sensitive to gas 
density n and the incident far-ultraviolet (FUV) flux Go 
( Stacey et aL"1|1991| ). Figure 2 shows ALESS65.1 compared 
with other low and high redshift galaxies, and solar metal- 
lically PDR model curves ( jKaufman et~al]p99| ). This dia- 
gram can be used to roughly estimate both n and Go for 
a galaxy, but there are some caveats, as outlined in |De| 
Breuck et aL] (|201l[>. These are: (i) the 12 CO luminosity 



is that of the ground state rotational line, (ii) [CII] emis- 
sion from the diffuse ionised medium and cosmic-ray-heated 
gas is assumed to be small compared to that from PDRs, 
and (iii) AGN and their related X-Ray Dissociation Regions 
(XDRs) are assumed to not contribute significantly to the 
FIR and [CII] luminosity. To be consistent with both De 
Breuck et al] ( |2011| ) and |Stacey et al.| ( [20ro| in Figure 2 we 
assume £co(2-i)/£co(i-o) = 7.2, which is 90% of its value 
if the gas was fully thermalised and optically thick. Cos- 
mic ray rates are greater in starbursts compared to normal 
galaxies but this does not seem to result in higher [CII] /CO 
ratios, so cosmic ray ionization does not appear to domi- 



nate the [CII] /CO ratio ( |De Breuck et al.||20iT| >, at least 
for local galaxies. ALESS65.1 is not detected in the 250ks 



Chan dra X-Ray observations of this region (Lehmer et al. 
|2005| , so it is not a luminous QSO (L 3 _44k e v < 2-3 xlO 44 
erg s~\ for jV H = - 10 23 5 cm" 2 ). The maximum likely 
AGN contribution to the FIR luminosity is estimated from 
a multicomponent (AGN and starburst) fit to the FIR us- 
ing Decompir software ( Mullaney et al.|2011 l and following 
the method described by ( Seymour et al.|2012 |. We use the 
24 /im flux density upper limit to constrain the short wave- 
length part of the IR SED and find the AGN contribution 
to the total FIR luminosity is <10%. ALESS65.1 therefore 
appears to be dominated by star-formation processes and 
the AGN contribution to [CII] and Lfir is minimal. 

In examining the PDR physical conditions, we multiply 



the 



CO(2-l) flux by a factor of two to account for the 

12/ 



line being optically thick, and hence only the CO emission 
coming from the illuminated PDR side is seen ( |Kaufman] 
et al.|1999||Hailey-Dunsheath et al.|2010[ ). We also multiply 



the [CII] flux by a factor of 0.7 to remove non-PDR con- 
tributions (e.g. |Hailey-Dunsheath et al.|[2010| |Stacey et al. 



20101. The CO geometry correction applies to all galaxies 
in Figure 2, and so does not affect the relative position of 
ALESS65.1 on the diagram compared to other galaxies. Us- 



ing the Kaufman et al. |l999j) models, we find ALESS65.1 
has Go ~ 10 d and n < lF~cm" 3 (Figure 2). Such a FUV 
radiation field is on the high end of those seen in low red- 
shift normal galaxies, but it is consistent with the strong 
FUV fields seen in local starbursts, nearby ULIRGs, and 
some z > 1 galaxies. This limit on Go and n implies a PDR 



temperature >300 K (Kaufman et al. 19991. Using Equa- 
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Table 1. Observed and derived properties of ALESS65.1 



Parameter 


Value 


Reference 


2[CII] 


4.4445 ± 0.0005 


S12 


/ [CH] 


5.4 ± 0.7 Jy km s" 1 


S12 


FWHM [CII] 


470 ± 35 km s" 1 


S12 


L [Cll\ 


(3.2 ± 0.4) x 10 9 L Q 


S12 


^FIR 


(1.38 ± 0.28) x 10 12 Lq 


S12 


^00(2-1) 


< 0.11 Jy km s- 1 (3a) 


this paper 


Af gas 


< 1.7 x 10 10 M (3a) 


this paper 


-^CO(2-l) 


< 8.5 X 10 6 L (3a) 


this paper 


T l 

^00(2-1) 


< 2.2 X 10 10 K km s' 1 pc 2 


this paper 



a converted from Lir(8 — 1000/xm) using Lpjjs (42 — 122 fim) 
= Ltr/1.45 (Stacey et al. 2010; De Breuck et al. 2011) 



tion 1 from Hailey-Dunsheath et al. (20101, we estimate the 



atomic gas associated with the PDR to be greater than 3 
X 10 9 Mq. We note that this atomic mass estimate is very 
uncertain, and it is highly dependent on n. The atomic gas 
mass M a ~ 2.9 x 10 9 + (3.2 X 10 12 /n), so for n = 10 3 cm" 3 
the atomic mass would be ~6 X 10 9 Mq. 

ALESS65.1 and the other two z > 4 SMGs shown 
in Figure 2, LESS J033229.4 (De Breuck et al. 2011) and 
HDF850.1 (Walter et al. 2012), have a similar FUV radi- 
ation field, Go, to local starbursts, but have much higher 
FIR luminosity, leading to suggestions that they are scaled- 
up versions of local starbursts. For a given Lfir the size 
of the emission region will increase for smaller Go- Fol- 
lowing Stacey et al. (2010), we scale up from M82 using 
two laws from Wolfire et al.l d 19901) to constrain the size: 



Go oc ALfir/D 3 if the mean free path of a UV photon A is 
small and Go oc Lyir/D 2 if the mean free path of a UV pho- 
ton is large. Applying these relations and using Go = 10 3 for 
ALESS65.1 yields a diameter of 1.1 - 2.1 kpc. This is con- 
sistent with the marginally resolved [CII] data which shows 
ALESS65.1 has a possible extent of 3.3 ± 1.7 kpc (S12). The 
same scaling relation applied to HDF850.1 results in a diam- 
eter of 1.8 - 4.6 kpc, which is consistent with the observed 
5.7 ± 1.9 kpc extent of the [CII] emission region (Walter et 
al. 2012). Similarly, LESS J033229.4 has an extent of ~ 4 
kpc (Coppin et al. 2010; De Breuck et al. 2011). Thus the 
starburst in all three z > 4 SMGs appears to be extended 
over galactic scales. 

Local starbursts and Galactic OB star forming regions 
lie on a line with [CII] /CO luminosity ratios of about 4100 
in Figure 2. Higher [CII]/CO ratios can also be found in 
low metallicity systems, such as 30 Doradus in the LMC, 
where the size of the [CII] emitting envelope of the cloud 
(relative to the the CO emitting core) is much larger than 
in more metal-rich systems ( Stacey et al.|1991 |. Metallicity 
may also affect the [CII] /CO ratios of the highest redshift 
galaxies. For example, LESS J033229.4 at z — 4.76 has a 
very high [CII] /CO ratio of « 10 4 and was initially thought 



to have sub-solar metallicity (De Breuck et al. 20111, but 



more recent [Nil] observations suggest that the metallicity 
of this galaxy is solar ( Nagao et al.|[2012 l. ALESS65.1 has 
L[cn]/ico > 2700, which is consistent with low metallic- 
ity, but it has an extent of several kpc so therefore the en- 
hanced [CII] emission is more likely to be due to the ex- 



tended star forming region. An extended (or less-dense) in- 
terstellar medium would have an increased fraction of UV 
photons available to ionize and excite gas, increasing the 



relative intensity of fine structure lines (e.g. Gracia-Carpio 
et al.|2011| | 



We have assumed a single-phase ISM in this work. 
Multi-phase ISMs are commonly required to ex plain the 
12 CO line ratios observed in local galaxies (e.g. Guesten 
et al. 1993 Ward et al. 2003[ ). In some cases a single- 



component ISM is unable to explain the 12 CO-excitation 
ladder (or spectral line energy distribution) observed in high 
redshift SMGs ([Carilli et al.|2010||Harris et al.|20Tol|Daniel- 



son et aL]|2011[ ). These studies found that the ISM in such 
sources is best described by a warm compact component sur- 
rounded by a cooler more extended one. Support for such a 
multi-component geometry in high-redshift sources has re- 
cently been found in the spatially-resolved 12 CO study of the 
gravitationally-magnified submm galaxy SMM J2135— 0102 
1 Swinbank et al.|2011 Danielson et al.|2011 1. Future higher- 
resolution continuum and line studies with ALMA may un- 
cover similar evidence in ALESS65.1, but the current data 
provide no direct evidence. 



for 20 hours to search for 
4.44 submillimetre galaxy. 



4 CONCLUSION 

We have observed ALESS65.1 
12 CO(2-l) emission in this z = 
We detect no 12 CO(2-l) emission in a spectrum reaching a 
rms sensitivity of 0.08 mjy beam -1 per 200 km s _1 channel. 

Adopting the FWHM from the ALMA detection of [CII] 
in ALESS65.1, we estimate a 3a limit to the 12 CO(2-l) lu- 
minosity of Lco(2-i) < 8.5 x 10 6 Lq and a cold gas mass 
upper limit of M gaa < 1.7 x 10 10 Mq. This implies a gas 
depletion timescale in ALESS65.1 of < 50 Myr, compara- 
ble to other z > 4 SMGs and consistent with this high red- 
shift population being the progenitors of z > 2 red-and-dead 
galaxies. 

We examine the physical conditions of the gas in 
ALESS65.1 using the L[cii]/iFm versus I/co(i-o)/£fir di- 
agram. We find ALESS65.1 has a strong FUV field compa- 
rable to local starbursts. The observed [CII] to FIR ratio, 
L[ch]/£fir = 2.3 x 10 , is high compared to local ULIRGs 
(as noted by S12). Combined with I/[cn]/£co > 2700, this 
high [CII] to FIR ratio is consistent with ALESS65.1 having 
more extended regions of intense star-formation than local 
ULIRGs. A possible, but less likely scenario, is ALESS65.1 
has low metallicity gas. 

Measurements of [CII] and 12 CO of a larger sample are 
needed to confirm whether z > 4 starbursts have enhanced 
[CII] emission compared to local galaxies, and whether this 
is because of metallicity effects, the relative size of PDR 
regions, a combination of the two, or other effects. Future 
surveys by ALMA will shed further light on the physical 
conditions of the gas in star forming galaxies in the early 
universe. 
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Figure 2. £[cii]/^fir versus Lco(l— o)/^fir f° r ALESS65.1 (red point with upper limit on £co(i— o)/-^fir) compared with Galactic 
star forming regions, starburst nuclei, normal galaxies, local ULIRGs, and high redshift (z > 1) sources. Empty symbols indicate low 
metallicity sources, which lie at high £fcii]/^'CO(l— 0)- Black lines represent the solar metallicity PDR model calculations for gas density 
(n) and FUV field strength (Go) from Kaufman et al. (1999). This figure is adapted from Stacey et al. (2010) with additional data from 
De Breuck et al. (2011) and Walter et al. (2012). ALESS65.1 has a higher Lr CI i]/Lpij^ ratio than that found in local ULIRGs, similar 
to the other two z > 4 SMGs, but its £[011] /^CO ratio is similar to local starbursts and other z > 1 sources. 
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